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Cavity-based large scale quantum information processing (QIP) may involve multiple cavities 
and require performing various quantum logic operations on qubits distributed in different 
cavities. Geometric-pbase-based quantum computing bas drawn much attention recently, 
which offers advantages against inaccuracies and local fluctuations. In addition, multiqubit 
gates are particularly appealing and play important roles in QIP. We here present a simple 
and efficient scheme for realizing a multi-target-qubit unconventional geometric phase gate 
in a multi-cavity system. This multiqubit phase gate has a common control qubit but differ¬ 
ent target qubits distributed in different cavities, which can be achieved using a single-step 
operation. The gate operation time is independent of the number of qubits and only two levels 
for each qubit are needed. This multiqubit gate is generic, e.g., by performing single-qubit 
operations, it can be converted into two types of significant multi-target-qubit phase gates 
useful in QIP and quantum Fourier transform. The proposal is quite general, which can be 
used to accomplish the same task for a general type of qubits such as atoms, NV centers, 
quantum dots, and qubits. 

M Ultiqubit gates are partieularly appealing and have been eonsidered as an attraetive build- 
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ing block for quantum information processing (QIP). There are two types of signifi¬ 
cant multiqubit gates, i.e., multiqubit gates with multiple control qubits acting on a single target 
qubit [1-6], and multiqubit gates with a single qubit simultaneously controlling multiple target 
qubits [7]. These two kinds of multiqubit gate play important roles in QIP. For instance, a multi- 
control-qubit gate has applications in quantum algorithms [8,9], error correction [10,11], and quan¬ 
tum Fourier transform [12]; and a multi-target-qubit gate is useful in quantum cloning [13], quan¬ 
tum algorithms [8,12,14], entanglement preparation [15], and error correction [16]. 

A multiqubit gate can in principle be constructed by using single-qubit and two-qubit ba¬ 
sic gates. However, when using the conventional gate-decomposition protocols to construct a 
multiqubit gate [17,18], the number of basic gates increases and the procedure usually becomes 
complicated as the number of qubits increases. Hence, building a multiqubit gate may become 
very difficult since each basic gate requires turning on and off a given Hamiltonian for a certain 
period of time, and each additional basic gate adds experimental complications and the possibility 
of more errors. During the past years, there is much interest in directly implementing multiqubit 
gates. Proposals have been presented for directly realizing both multi-control-qubit gates [1-6] and 
multi-target-qubit gates [7] in various physical systems. Note that the gate implementation using 
these previous proposals [1-7] was based on non-geometric dynamical evolution. 

In recent years, much attention has been paid to fault tolerant quantum computing, which 
operates essentially based on geometric phase [19-22]. The geometric phase is determined by the 
global features of the evolution path, which offers potential advantages against local fluctuations. 
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A great deal of work has been devoted to construct geometric phase gates, which are usually 
divided into two categories: conventional geometric phase gates and unconventional geometric 
phase gates. The former needs to remove the dynamical phase by choosing the cyclic evolution in 
dark states [23] or employing multi-loop schemes (the evolution is driven by a Hamiltonian along 
several closed loops) [24,25]. But the latter does not require performing additional operations 
to cancel the dynamical phase, because the total phase is dependent only on global geometric 
features and independent of initial states of the system [26,27]. The unconventional geometric 
phase gate was first proposed in [26]. According to [26], an unconventional geometric phase gate 
is characterized by a unitary operator [/({y}), where 7 is the total phase, which consists of a 
geometric phase and a dynamic phase (see [26]). 

During the past years, proposals for realizing both conventional geometric phase gates [23- 
25,28-31] and unconventional geometric phase gates [26,27,32] were presented. Moreover, high fi¬ 
delity conventional geometric phase gates have been experimentally demonstrated using NMR [33], 
superconducting qubits [34], and nitrogen-vacancy (NV) centers [35,36]. In addition, a two-qubit 
unconventional geometric phase gate has been experimentally realized in a trapped ion system [37]. 
However, previous works are mainly focused on how to construct single-qubit or two-qubit con¬ 
ventional/unconventional geometric phase gates [23-37], or implementing a multi-control-qubit 
gate [1-6] and a multi-target-qubit gate [7] based on non-geometric dynamical evolution. 

In this work, we consider how to implement a multi-target-qubit unconventional geometric 
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phase gate, which is described by the following transformation: 

n n 

i=i i=i 

n n 

i=i i=i 

where subscript A represents a control qubit, subscripts (1, 2, n) represent n target qubits (1, 2, n), 

n 

and n Kj) (with ij e {+, —}) is the n-target-qubit computational basis state. For n target qubits, 

i=i 

there are a total number of 2” computational basis states, which form a set of complete orthogonal 
bases in a 2"^-dimensional Hilbert space of the n qubits. Equation (1) shows that when the con¬ 
trol qubit A is in the state |-|-) (|—)), a phase shift happens to the state |+) (|—)) but nothing 
happens to the state |—) (|-|-)) of the target qubit j (j = 1, 2,..., n). For instance, under the trans¬ 
formation (1), one has: (i) the state transformation described by following Eq. (18) for a two-qubit 
phase gate on control qubit A and target qubit j, and (ii) the state transformation described by 
Eq. (21) below for a three-qubit phase gate on control qubit A and two target qubits (1, 2). Note 
that the multiqubit phase gate described by Eq. (1) is equivalent to such n two-qubit phase gates, 
i.e., each of them has a common control qubit A but a different target qubit 1, 2,..., or n, and 
the two-qubit phase gate acting on the control qubit A and the target qubit j (j = 1, 2,..., n) is 
described by Eq. (18). 

The multiqubit gate described by Eq. (1) is generic. For example, by performing a single- 

n 

qubit operation such that |-|-)^ ^ H and |—^ )j but nothing to |—)a and 

i=i 

the transformation (1) becomes 

n n 

I+)a]^K, 7) |+)A]^Ki), 

i=i j=i 
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( 2 ) 


i=i i=i 

which implies that when and only when the control qubit A is in the state |—), a phase shift 
happens to the state |—) of the target qubit j but nothing otherwise (Fig. 2). For dj = n 12^ the state 
transformation (2) corresponds to a multi-target-qubit phase gate, i.e., if and only if the control 
qubit A is in the state |—), a phase flip from the sign -f to — occurs to the state |—) of each 
target qubit. Such a multiqubit phase gate has many applications in QIP [13-16]. In addition, for 
= n/2\ the state transformation (2) corresponds to a multi-target-qubit phase gate, i.e., if and 
only if the control qubit A is in the state |—), a phase shift dj = 7r/2-^ happens to the state |—) of 
each target qubit. It is noted that this multi-target-qubit gate is equalivent to a multiqubit gate with 
different control qubits acting on the same target qubit (Fig. 3), which is a key element in quantum 
Fourier transform [8,12]. 

In what follows, our goal is propose a simple method for implementing a generic uncon¬ 
ventional geometric (UG) multi-target-qubit gate described by Eq. (1), with one qubit (qubit A) 
simultaneously controlling n target qubits (1, 2,..., n) distributed in n cavities (1, 2,..., n). We be¬ 
lieve that this work is also of interest from the following point of view. Large-scale QIP usually 
involves a number of qubits. Placing many qubits in a single cavity may cause some fundamental 
problems such as introducing the unwanted qubit-qubit interaction, increasing the cavity decay, 
and decreasing the qubit-cavity coupling strength. In this sense, large-scale QIP may need to place 
qubits in multiple cavities and thus require performing various quantum logic operations on qubits 
distributed in different cavities. Hence, it is important and imperative to explore how to realize 
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multiqubit gates performed on qubits that are spatially-separated and distributed in different cavi¬ 


ties. 


As shown below, this proposal has the following features and advantages: (i) The gate oper¬ 
ation time is independent of the number of qubits; (ii) The proposed multi-target-qubit UG phase 
gate can be implemented using a single-step operation; (iii) Only two levels are needed for each 
qubit, i.e., no auxiliary levels are used for the state coherent manipulation; (iv) The proposal is 
quite general and can be applied to accomplish the same task with a general types of qubits such as 
atoms, superconducting qubits, quantum dots, and NV centers. To the best of our knowledge, this 
proposal is the first one to demonstrate that a multi-target-qubit UG phase gate described by (1) 
can be achieved with one qubit simultaneously controlling n target qubits distributed in n cavities. 

In this work we will also discuss possible experimental implementation of our proposal and 
numerically calculate the operational fidelity for a three-qubit gate, by using a setup of two super¬ 
conducting transmission line resonators each hosting a transmon qubit and coupled to a coupler 
transmon qubit. Our numerical simulation shows that highly-fidelity implementation of a three- 
qubit (i.e., two-target-qubit) UG phase gate by using this proposal is feasible with present-day 
circuit QED technique. 

Results 

Model and Hamiltonian. Consider a system consisting of n cavities each hosting a qubit and 
coupled to a common qubit A [Fig. 1(a)]. The coupling and decoupling of each qubit from its 
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cavity can be achieved by prior adjustment of the qubit level spacings. For instance, the level 
spacings of superconducting qubits can be rapidly adjusted by varying external control parameters 
(e.g., magnetic flux applied to the superconducting loop of a superconducting phase, transmon, 
Xmon or flux qubit; see, e.g., [38-41]); the level spacings of NV centers can be readily adjusted by 
changing the external magnetic field applied along the crystalline axis of each NV center [42,43]; 
and the level spacings of atoms/quantum dots can be adjusted by changing the voltage on the 
electrodes around each atom/quantum dot [44]. The two levels of coupler qubit A are denoted 
as \g)A and \e)A while those of intracavity qubit j as \g)j and \e)j (j = 1, 2, • • •, n). Applying 
a classical pulse to qubit A and a classical pulse to each intracavity qubit j [Fig. l(b),(c)]. For 
identical qubits, we have cc = UegA = where a; is the pulse frequency and (coegj) is the 
\g) -H- |e) transition frequency of qubit A (qubit j). The system Hamiltonian in the interaction 
picture reads (in units of h= 1) 

n n 

Hi = '^ gj{e~''^^^a]a~ + h.c.) + ^ gA^ + h.c.) 

i=i i=i 

n 

+ ^ ( 3 ) 

i=i 

where a] is the photon creation operator for the mode of cavity j, = |e)A(fi'| (\e)j{g\) 

and <j^ (a~) = |5')A(e| (\g)j{e\) are the raising and lowering operators for qubit A (qubit j), 
Sj = Uegj — oJcj and 5^^ = ooeg^ — are detunings (with Uc^ being the frequency of cavity j), H 
is the Rabi frequency of the pulse applied to each qubit, gAj (gj) is the coupling constant of qubit 
A (j) with cavity j. We choose |±)j = {\e)j ± \g)j)/\/2 and |±)a = (|e)A ± |5 ')a)/\/2 as the 
rotated basis states of qubit j and qubit A, respectively. 
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In a rotated basis {I+) one has + o'; ) /2and(T; = — c^i) /2, 

where = |+)i(+l “ |-)z(-|’3"+ = |+)/(-|, and a" = |-)K+I- Here, / = 1,2,3, •••n, A 
Hence, the Hamiltonian (3) can be expressed as 


1=1 
” 1 

+ - ^ a ) + h.c] 

n 

+ . (4) 

1=1 


n 


In a new interaction picture under the Hamiltonian Hq= Yh , one obtains from Eq. (4) 


1=1 




1=1 




2^Aje - e + h.c]. 


( 5 ) 


1=1 

In the strong driving regime 2f2 3> {gj, |(5j | jfi'Aj, I'^aJ}, one can apply a rotating-wave approx¬ 
imation and eliminate the terms that oscillate with high frequencies. Thus, the Hamiltonian (5) 
becomes 


H'j = 


n 

E 

j=i 


ri/lA 


+ h.c.) + 


n 


1 


7. 

1=1 


—i&At t 

.e J a' 


+ h.c.). 


For simplicity, we set 


( 6 ) 


9Aj = gj, 5j = 6 aj- (7) 

The first term of condition (7) can be achieved by adjusting the position of qubit j in cavity j, and 
second term can be met for identical qubits. Thus, the Hamiltonian (6) changes to 

n 

1=1 
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with 


= ^9j{e + a,^), (9) 

where Heff,j is the effective Hamiltonian of a subsystem, which consists of qubit A, intracavity 
qubit j, and cavity j. In the next section, we first show how to use the Hamiltonian (9) to construct 
a two-qubit UG phase gate with qubit A controlling the target qubit j. We then discuss how to use 
the effective Hamiltonian (8) to construct a multi-qubit UG phase gate with qubit A simultaneously 
controlling n target qubits distributed in n cavities. 

Implementing multiqubit UG phase gates. Consider a system consisting of the coupler 
qubit A and an intracavity qubit j, for which |±)j (|±)a) are eigenstates of the operator (o'za) 
with eigenvalues ±1. In the rotated basis {|+)A|+)i, |+)a|—) i, |“)A|+)i, |“)a|—) i}, the Hamil¬ 
tonian (9) can be rewritten as 

= 9j{e x ^|-f)A|+)j (+Ia (+lj “ l“)A|-)i (-|a (~\j) > (1*^) 

and thus the time evolution operator UAj{t) corresponding to the Hamiltonian can be ex¬ 

pressed as 

UAj{t) = ^+-r,j(f)l+)A|+)i (+Ia (+lj + l+)A|-)i (+Ia (“Ij 

+ |-)A|+)i (-|a (+lj + ^—,i(^)l“)A|-)i (-|a (“Ij ) (11) 

where {t) and U _ (t) are given by 

Upp,j{t) = fjexp{-i / Hpp^^{T)dT) 

Jo 
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= Tjexp[—igjepp [ + e^^^'^aj)dT] 

Jo 

N 

= lim exp[—+ e*^^''^’*aj)Ar] 

n=l 

N 

n=l 

= D(/ciapp,,)e*^--, (12) 


with 

Hpp,, {t) = {p\^ {p\. Heff,^ \p)A\p)j = gjEppie-^^^^a] + (13) 

where pp G {++, —}, p G {+, —}, e++ = —e _= 1, D is the displacement operator (for de¬ 

tails, see Methods below), Tj is the time ordering operator and At = t/N is the time interval. From 
Eq. (12) and Eq. (31) below, one obtains dappj = —igjeppe~'^^^'^dT and Opp^ = Irn{J^ app jdappj). 
Thus, one has 


appj = I dappj = - 1 ), 


^pp,j 




^ppd = —T / 0-— cos 5jt)dt, 

Oj Jo 


(14) 


where Tj is the evolution time required to complete a closed path. 


If t = Tj is equal to 2mjT/\5j\ with a positive integer rrij, we have j\,oipp,j = 0 according 
to Eq. (14), which shows that when cavity j is initially in the vacuum state, then cavity j returns 
to its initial vacuum state after the time evolution completing a closed path. Thus, it follows from 
Eq. (12) that we have 

Upp,^ (Tj ) = D (0) . (15) 
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Here Opp^ is the total phase given by Eq. (14), which is acquired during the time evolution from 
t = 0 to t = Tj. Note that 9ppj consists of a geometric phase and a dynamical phase. 


It follows from Eqs. (11) and (15) that the cyclic evolution is described by 

UAj{Tj) = e*®++’-’|+)A|+)j (+U (+|j + |+)A|-)i (+U (-|j 

+ |-)A|+)i (-U (+lj + (-U (“Ij • (16) 


Eq. (14) shows that 6ppj is independent of index pp. Thus, we have = 6 _j = 6j. Eurther, 

according to Eq. (14), after an integration for Tj = 2mj7r/\6j \ (set above), we have 


e, = -fr, = 


2mjTTg^. 


which can be adjusted by varying the coupling strength gj and detuning 6j. Note that a negative 
detuning 6j < 0 (see Eig. 1) has applied to the last equality of Eq. (17). The unitary operator (16) 
describes a two-qubit UG phase gate operation. Eor 6j ^ 2mT with an integer n, the phase gate is 
nontrivial. After returning to the original interaction picture by performing a unitary transformation 
U = exp{—-f VtazA)Tj}, we obtain the following state transformations: 1+)^!+)^ ^ 

^ ^ |-)^|+)„ and |-)^|-), ^ 

|—)^|—)j, which can be further written as 


|+)A|+)i ^ e*®-’l+)A|+)i 
l+)A|-)i ^ l+)A|-)i 

l“)A|-)i (18) 
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where we have set QTj = kn (k is a positive integer). For Tj = 2mjn/\Sj\, we have 2Q = 
k\6j\/mj. The result (18) shows that a two-qubit UG phase gate was aehieved after a single-step 
operation deseribed above. 

Now we expand the above proeedure to a multiqubit ease. Consider qubit A and n qubits (1,2, 
• • •, n) distributed in n eavities [Fig. 1(a)]. From Eqs. ( 8 ) and (9), one ean see that: (i) eaeh term 
of Heff aets on a different intra-eavity qubit but the same eoupler qubit A, and (ii) any two terms 
of Heff, eorresponding to different j, eommute with eaeh other: [iFe//,j,= 0 (j 7 ^ /c = 
1, 2, • • • n). Thus, it is straightforward to show that the eyelie evolution of the eavity-qubit system 
is deseribed by the following unitary operator 

n 

U(T) = 1[UaAT,), (19) 

J = 1 

where UAj{Tj) is the unitary operator given in Eq. (16), whieh eharaeterizes the eyelie evolution 
of a two-qubit subsystem (i.e., qubit A and intraeavity qubit j) in the rotated basis |-l-)yi|+)j, 
l+)A|-)j, |-)A|+)i, and |-)yi|-)j. 

By ehanging the detunings 5j (e.g., via prior design of eavity j with an appropriate fre- 
queney), one ean have 

mi/ 6 i = 7712/62 =,■■■, = ^n/Sn, (20) 

whieh leads to Ti = T 2 =,•••,= = T, i.e., the evolution time for eaeh of qubit pairs {A, 1), 

{A, 2), ■ • •, and {A,n) to eomplete a eyelie evolution is identieal. Eor the setting here, we have 
9j = —-^T resulting from Eq. (17). Henee, one ean easily find from Eqs. (18) and (19) that 
after a eommon evolution time T, the n two-qubit UG phase gates eharaeterized by a jointed 


12 



unitary operator U(T) of Eq. (19), which have a common control qubit A but different target qubits 
(1, 2, n), are simultaneously implemented. As discussed in the introduction, the n two-qubit UG 
phase gates here are equivalent to a multiqubit UG phase gate described by Eq. (1). Hence, after the 
above operation, the proposed multiqubit UG phase gate is realized with coupler qubit A (control 
qubit) simultaneously controlling n target qubits ( 1 , 2 , • • •, n) distributed in n cavities. 

To see the above more clearly, consider implementing a three-qubit (two-target-qubit) UG 
phase gate. Eor three qubits, there are a total number of eight computational basis states, denoted 
by {|+)a|+)i|+) 2 , |+)a|+)i|-) 2 , •••, |-)a|-)i|-) 2 }- According to Eqs. (18) and (19), one 
can obtain a three-qubit UG phase gate, which is described by 


+ )a +)l +)2 

e*(^^+'^)|+)^|+)i|+)2, 

-)a +)l +)2 ^ -)a +)l +)2, 

+ )a +)l -)2 ^ 

+)a +)i -)2, 

-)a +)l -)2 -)a +)l -)2, 

+ )a -)l +)2 

+)a -)i +)2, 

-)a -)i +)2 ^ -)a -)i +)2, 

+ )a -)l “)2 ^ 

+)a -)i “)2, 

|-)A|-)i|-)2^e*(^^+^^)|-)^|-)i|-)2. (21) 


As discussed in the introduction, by applying single-qubit operations, this three-qubit UG phase 
gate described by Eq. (21) can be converted into a three-qubit phase gate illustrated in Eig. 2 or 
Eig. 3 for n = 2. In the next section, as an example, we will give a discussion on the experimental 
implementation of this three-qubit UG phase gate for the case of 6 *i = 6*2 = tt/2. Based on Eq. 
(17) and for Ti = T 2 (see above), one can see that the 6*1 = 62 corresponds to gf/Si = (y'|/( 52 , 
which can be met by adjusting gj (e.g., varying the position of qubit j in cavity j) or detuning 5j 
(e.g., prior adjustment of the frequency of cavity j) (j = 1 , 2 ). 
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Possible experimental implementation. Superconducting qubits are important in QIP due 
to their ready fabrication, controllability, and potential scalability [38,45,46]. Circuit QED is ana¬ 
logue of cavity QED with solid-state devices coupled to a microwave cavity on a chip and is 
considered as one of the most promising candidates for QIP [45-49]. In above, a general type of 
qubit, for both of the intracavity qubits and the coupler qubit, is considered. As an example of 
experimental implementation, let us now consider each qubit as a superconducting transmon qubit 
and each cavity as a one-dimensional transmission line resonator (TER). We consider a setup in 
Eig. 4 for achieving a three-qubit UG phase gate. To be more realistic, we consider a third higher 
level I/) of each transmon qubit during the entire operation because this level |/) may be excited 
due to the |e) -H- |/) transition induced by the cavity mode(s), which will affect the operation 
fidelity. Erom now on, each qubit is renamed “qutrit” since the three levels are considered. 

When the intercavity crosstalk coupling and the unwanted |e) -H- |/) transition of each qutrit 
are considered, the Hamiltonian (3) is modified as follows 

hi = Hi + Qi, (22) 

where Hj is the needed interaction Hamiltonian in Eq. (3) for n = 2, while 0/ is the unwanted 

interaction Hamiltonian, given by 

2 _ 2 __ 

0/ = + h.c^ + '^9Aj + h.c^ + gi 2 (e*^*aia^ -f h.c.) 

j=i i=i 

2 

+ ^D , + h.c.] + n ^ + h.c.], (23) 

i=i 

where = \f)j (e| and <7^^^ = |/)^ (e| . The first term describes the unwanted off-resonant 
coupling between cavity j and the |e) -H- |/) transition of qutrit j, with coupling constant gj 
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and detuning 5j = ujfej — <^cj [Fig- 5(a,b)], while the second term is the unwanted off-resonant 
coupling between cavity j and the |e) -H- |/) transition of qutrit A, with coupling constant gAj and 
detuning Saj = Ufe^^—Ucj [Fig. 5(c)]. The third term of Eq. (23) describes the intercavity crosstalk 
between the two cavities, where A = Uc 2 — cua =61 — 82 is the detuning between the two-cavity 
frequencies and gi 2 is the intercavity coupling strength between the two cavities. The last two 
terms of Eq. (23) describe unwanted off-resonant couplings between the pulse and the |e) -H- |/) 
transition of each qutrit, where is the pulse Rabi frequency. Note that the Hamiltonian (23) 
does not involves \g) -H- |/) transition of each qutrit, since this transition is negligible because of 
cUc., CJ < Ufg.,Ufg^ O' = 1,2) (Eig. 5). 

When the dissipation and dephasing are included, the dynamics of the lossy system is deter¬ 
mined by the following master equation 

i=i 

+ [^Jei] + ^fgi^ [^/sJ } 

1=1,2,A 

1=1,2,A 

T ^ ^ {F^^i^e (^eeiP^ee; ^eeiP/2 P^eei/2)} , (24) 

l=l,2,A 

where = \g)^ {f\, = \e)i (e| , a//, = |/)^ (/|; and £ [A] = ApA+ - A+Ap/2 - pA+A/2, 

with A = ttj, af, '^fei^'^Jgi- Here, kj is the photon decay rate of cavity aj 0 = 1, 2). In addition, 
F; is the energy relaxation rate of the level |e) of qutrit /, F/g, (F fg^) is the energy relaxation rate of 
the level |/) of qutrit / for the decay path |/) ^ |e) (|p)), and F; (rz,(p/) is the dephasing rate of 
the level |e) (|/)) of qutrit / (/ = 1, 2, A). 
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The fidelity of the operation is given by 


^ ( 25 ) 

where \ipid) is the output state of an ideal system (i.e., without dissipation, dephasing, and crosstalk 
considered), while p is the final density operator of the system when the operation is performed in 
a realistic physical system. As an example, we consider that qutrit I is initially in a superposition 
state l/\/2 (|+); + |—);) (/ = 1, 2, A) and cavity 1 (2) is initially in the vacuum state. In this case, 
we have = Wid) ® | 0 )ci| 0 )c 2 , where 

= (1/\/8)(—|+)a|+)i|+)2 + *|+)a|+)i|~)2 + *|+)a|~)i|+)2 + |+)a|~)i|~)2 
+ |~)a|+)i|+)2 + *|~)a|+)i|~)2 +*|~)a|~)i|+)2 ~ |~)a|~)i|~)2)5 (26) 

which is obtained based on Eq. (21) and for 6 *i = 6*2 = 7 r/ 2 . 

We now numerically calculate the fidelity of the gate operation. Without loss of generality, 
consider identical transmon qutrits and cavities. Setting mi = 1 and m 2 = 2, we have 82 = 
25i because of Eq. (20), which corresponds to gi /§2 = l/\/2 for 6*1 = 62 - In order to satisfy 
the relation 2f2 3> 1^21 and 2VL = A;|( 52 |/ 2 , we set A; = 12. In addition, we have pj 
gAj ~ V^gAj {j = 1,2), and (l ~ \/2f2 for the transmon qutrits [50]. Eor a transmon qutrit, 
a ratio 5% of the anharmonicity between the \g) -H- |e) transition frequency and the |e) -H- |/) 
transition frequency is readily achieved in experiments. Thus, we set 5j = 5j — O.ObUegj and 
^Aj = 8 j — O.ObUeg^ (j = 1,2). Eor transmon qutrits, the typical transition frequency between 
two neighbor levels is between 4 and 10 GHz [51]. Therefore, we choose uj^gj2'K^ ujegj2'K ~ 6.5 
GHz. Other parameters used in the numerical calculation are as follows: = 10.0 ps. 


16 



^ = 30.0 /iS, = 11.5 /iS, = 45.0 US (/ = 1, 2, A), and Kj ^ = 15.0 ns (j = 1, 2). 

To test how the inter-cavity crosstalk affects the gate fidelity, we plot Fig. 6 for gi 2 = 
0, O.lgfi, 0 . 25 ( 1 , 0 . 351 , which shows the fidelity versus For simplicity, the dissipation and 

dephasing of the system are not considered in Fig. 6. As depicted in Fig. 6, the effect of the 
inter-cavity coupling is negligible as long as 512 < O.I 51 . 

Figure 7 shows the fidelity versus (5i/27r, which is plotted by setting 512 = O.I 51 and now 
taking the systematic dissipation and dephasing into account. From Fig. 7, one can see that for 
(5i/27r = —3.57 MHz, a high fidelity 97.0% is achievable for a three-qubit UG phase gate. For 
(5i/27r = -3.57 MHz, we have T = Ti = T 2 = 0.28/iS, 5i/27r = 1.79 MHz, and 52/271 = 2.52 
MHz. The values of 51 and 52 here are readily available in experiments [52]. 

The condition 512 < O.I51 is easy to satisfy with the cavity-qutrit capacitive coupling shown 
in Fig. 4. When the cavities are physically well separated, the inter-cavity crosstalk strength is 
512 ~ 5 AiC' 2 /Gs, 5 yi 2 C'i/C's, where Gs = Ci + C 2 + Cq {Cq is the qutrit’s self-capacitance) 
[53,54]. For Gi, G2 ~ 1 fF and Gs ~ 100 fF (typical values in experiments), one has 512 = O.OI51. 
Thus, the condition 512 < O.I51 is readily achievable in experiments. 

Energy relaxation time Ti and dephasing time T 2 of the level |e) can be made to be on 
the order of 20 — 60 /iS for state-of-the-art transom devices [55]. For transmon qutrits, we have 
the energy relaxation time T[ ~ Ti /2 and dephasing time T/ ~ T 2 of the level |/), which are 
comparable to Ti and T 2 , respectively. With Ueg^l2TT, ujegJ^Ti ~ 6.5 GHz chosen above, we have 
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oJci/^Ti ~ 6.503 GHz and uJc 2 /‘^t^ ~ 6.507 GHz. For the cavity frequencies here and the values 
of ^7^ and K 2 ^ used in the numerical calculation, the required quality factors for the two cavities 
are Qi ~ 6.126 x 10® and Q 2 ~ 6.130 x 10®. Note that superconducting coplanar waveguide 
resonators with a loaded quality factor Q ~ 10® were experimentally demonstrated [56,57]. We 
have numerically simulated a three-qubit circuit QED system, which shows that the high-fidelity 
implementation of a three-qubit UG phase gate is feasible with current circuit QED technique. 

Discussion 

A simple method has been presented to realize a generic unconventional geometric phase gate of 
one qubit simultaneously controlling n spatially-separated target qubits in circuit QED. As shown 
above, the gate operation time is independent of the number n of qubits. In addition, only a single 
step of operation is needed and it is unnecessary to employ three-level or four-level qubits and not 
required to eliminate the dynamical phase, therefore the operation is greatly simplified and the ex¬ 
perimental difficulty is significantly reduced. Our numerical simulation shows that highly-fidelity 
implementation of a two-target-qubit unconventional geometric phase gate by using this proposal 
is feasible with the present circuit QED technique. The proposed multiqubit gate is generic, which, 
for example, can be coonverted into two types of important multi-target-qubit phase gates useful 
in QIP and quantum Eourier transform. This proposal is quite general and can be applied to ac¬ 
complish the same task with various types of qubits such as atoms, quantum dots, superconducting 
qubits, and NV centers. 
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Methods 


Geometric phase. Geometric phase is induced due to a displacement operator along an arbitrary 
path in phase space [58,59]. The displacement operator is expressed as 

D{a) = (27) 

where and a are the creation and annihilation operators of an harmonic oscillator, respectively. 
The displacement operators satisfy 

D{ai)D{a2) = D{ai + (28) 


For a path consisting of N short straight sections Aa^, the total operator is 


Dt = D{Aajsf) ■ ■ ■ D{Aaj) 

N N j-1 

= DC^Aaj)exp[iImC^Aaj'^Aal)]. 

j=l j=2 k=l 


(29) 


An arbitrary path c can be approached in the limit N ^ oo. Therefore, Eq. (29) can be rewritten 


as 


Dt = D{ da)e 


,ie 


(30) 


with 


0 = Im{ / a*da). 


(31) 


For a closed path, we have 


Dt = D{0)D® = e*®, 


(32) 
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where 0 is the total phase whieh eonsists of a geometrie phase and a dynamieal phase [22]. In 
above, equations (27-32) have been adopted for realizing an UG phase gate of one qubit simulta¬ 
neously eontrolling n target qubits. 
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Figure 1: (a) Diagram of a coupler qubit A and n cavities each hosting a qubit. A blue 
square represents a cavity while a green dot labels a qubit placed in each cavity, which can be 
an atom or a solid-state qubit. The coupler qubit A can be an atom or a quantum dot, and can 
also be a superconducting qubit capacitively or inductively coupled to each cavity, (b) Cavity 
j is dispersively coupled to qubit j (placed in cavity j) with coupling constant gj and detuning 
Sj < 0. (c) The coupler qubit A dispersively interacts with cavity j, with coupling constant gAj 
and detuning Saj < 0 (j = 1, 2,..., n). Here, Saj = Sj, which holds for identical qubits A and j. 

Figure 2: (a) Schematic circuit of a phase gate with qubit A (a black dot) simultaneously 
controlling n target qubits (squares), (b) This multiqubit phase gate illustrated in (a) consists of 
n two-qubit phase gates, each having a shared control qubit (qubit A) but a different target qubit 
(qubit 1, 2, • • • ,orn). Here, the element 26j represents a phase shift exp(i26*j), which happens to 
the state |—) of target qubit j (j = 1, 2,..., n) when and only when the control qubit A is in the 
state |—) but nothing happens otherwise. For 29j = tt, this gate corresponds to a multi-target-qubit 
phase gate (usedful in QIP [13-16]), i.e., if and only if the control qubit A is in the state |—), a 
phase flip from the sign -f to — occurs to the state |—) of each target qubit. 

Figure 3: Schematic circuit of the n successive two-qubit phase gates in quantum Fourier 
transform. Here, each two-qubit phase gate has a shared target qubit (qubit A) but a different 
control qubit (qubit 1, 2, • • • ,orn). The element 71/2^ represents a phase shift exp(i7r/2'^), which 
happens to the state |—) of target qubit A if and only if the control qubit j is in the state |—) (j = 
1, 2,..., n). For any two-qubit controlled phase gate described by the transformation |-l-)A|+)j ^ 
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|+)A|+)i, I+)a|—) i ^ 1“)^!+)^, |“)A|+)j l“)A|+)i, and |—)a|—) i —^ it is 

clear that the roles of the two qubits can be interchanged. Namely, the first qubit can be either the 
control qubit or the target qubit, and the same applies to the second qubit. When the second (first) 
qubit is a control qubit, while the first (second) qubit is a target, the phase of the state |—) of the 
first (second) qubit is shifted by when the second (first) qubit is in the state |—), while nothing 
happens otherwise. Thus, the quantum circuit here is equivalent to the circuit illustrated in Fig. 2 
for20, =7r/2Mj = l,2,...,n). 

Figure 4: Setup of two cavities (1,2) connected by a superconducting transmon qubit A. 
Here, each cavity represents a one-dimensional coplanar waveguide transmission line resonator, 
qubit A is capacitively coupled to cavity j via a capacitance Cj (j = 1,2). The two green dots 
indicate the two transmon qubits (1,2) embedded in the two cavities, respectively. The interaction 
of qubits (1,2) with their cavities is illustrated in Figs. 5(a) and (b), respectively. The interaction of 
qubit A with the two cavities is shown in Fig. 5(c). Due to three levels for each qubit considered 
in our analysis, each qubit is renamed as a qutrit in Fig. 5. 

Figure 5: Schematic diagram of qutrit-cavity interaction, (a) Cavity 1 is coupled to the 
\g) -H- |e) transition with coupling strength gi and detuning 5i, but far-off resonant with the |e) -H- 
I/) transition of qutrit 1 with coupling strength gi and detuning (5i. (b) Cavity 2 is coupled to 
the \g) -H- |e) transition with coupling strength g 2 and detuning <52, but far-off resonant with the 
|e) -H- I/) transition of qutrit 2 with coupling strength g 2 and detuning S 2 . (c) Cavity 1 (2) is 
coupled to the l^f) -H- |e) transition of qutrit A with coupling strength gA^ (fl'Aa) and detuning 5ai 
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(5^2)’ but far-off resonant with the |e) -H- |/) transition of qutrit A with eoupling strength 
0 A 2 ) and detuning Sa, (Sa^)- Here, Sj = - Uc^,5j = ujej - ^c^^a^ = ujegA - and 

5Aj = ujfsA - ^Cj O' = 1, 2), where Ueg^ (ufej) is the \g) -H- |e) (|e) -H- |/)) transition frequeney 
of qutrit j, Ueg^ is the \g) <H- |e) (|e) -H- |/)) transition frequeney of qutrit A, and Uc^ is the 

frequeney of eavity j. 

Figure 6: Fidelity versus 5i/2ti^ plotted for different intereavity eoupling strengths but with¬ 
out eonsidering the systematie dissipation and dephasing for simplieity. 

Figure 7: Fidelity versus (5i/27r, plotted for gi 2 = 0 and gi 2 = O.lgfi and by taking the sys¬ 
tematie dissipation and dephasing into aeeount. The parameters used in the numerieal simulation 
for Figs. 6 and 7 are referred to the text. 
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